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Using  the  feres  and  moment  expressions,  s tuple  an  chan  leal  sys  tests  equivalent 


All  of  the  above  developments  assume  a non-vlseous  Incompressible 
fluid  subject  to  linearised  boundary  conditions.  All  tank  actions  (except 
that  nocraal  to  the  n»*m  free  surface)  are  restricted  to  snail  accelerations 


tion. 


Codg2. 

.d/or 


There  has  been  increasing  interest  recently  in  aircraft  dynamics 
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problems  which  are  affected  by  the  motion  of  fuel  in  partially  full  tanks  ’ . 
The  stability  and  response  of  a rigid  airplane  is  a problem  of  this  type. 

It  is  intuitively  apparent  that  a large  mass  of  fuel  with  a free  surface 
may  have  serious  dynamic  effects  on  the  airplane , particularly  if  there  is 
a possibility  of  resonance. 

Available  Information 

The  understanding  of  these  problems  requires  first  an  understanding 
of  the  dynamics  of  fuel  motion.  A great  deal  of  the  necessary  information 
has  been  available  for  some  time^*,  particularly  the  nature  of  the  free 
oscillations  of  liquid  in  some  tanks  of  simple  shape.  Extension  of  this 
material  to  include  forced  oscillations  makes  it  possible  to  determine  the 
forces  communicated  to  the  tank  walls  by  an  arbitrary  motion  of  the  tank 
(or  airplane).  Prior  work*  in  this  direction  initiated  by  Stewart  and 
Lorell  was  in  part  responsible  for  the  present  investigation.  Similar 
problems  have  been  investigated  in  the  field  of  seismology.  Ref.  5 gives 
experimental  results  and  additional  references  to  theoretical  and  experimental 
work  in  this  field. 

Methods  Used 

For  convenience  in  solving  the  complete  dynamics  problem  it  may  be 
desirable  to  replace  the  fuel  by  a simple  mechanical  system.  Therefore, 
in  addition  to  presenting  the  Laplace  transforms  of  forces  and  moments, 
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this  report  suggests  mechanical  systems  which  produce  the  same  forces  and 
moments  as  the  oscillating  fuel  when  the  container  is  given  an  arbitrary 
motion.  (The  "arbitrary"  motion  must,  of  course,  fall  within  the  scope  of 
the  linearization  processes  used.)  In  a previous  report^  the  representation 
of  fuel  by  an  equivalent  pendulum  was  discussed.  The  characteristics  of  the 
equivalent  pendulum  were  determined  only  from  consideration  of  fluid  motions 
in  a stationary  tank.  Ibis  was  adequate  to  define  the  frequency  and  the 
product  of  mass  times  amplitude  of  motion  for  the  pendulum,  so  that  forces 
corresponding  to  a given  amplitude  of  the  fundamental  wave  motion  were  dupli- 
cated. Ibis  provided  sufficient  information  for  the  study  of  Borne  stability 
problems  as  affected  by  fuel  "sloshing."  The  present  report  offers  additional 
information  so  that  the  more  general  problem  of  aircraft  response  to  arbitrary 
disturbances  can  be  considered,  including  the  effect  of  angular  tank  motions. 


PROCEDURE 


Assumptions 

In  order  to  simplify  the  problem  a number  of  restrictions  are  applied 
concerning  the  nature  of  the  fluid,  the  motions  of  the  fluid  and  the  shape 
of  the  tank.  [Hie  fluid  is  assumed  incompressible  (which  requires  no  dis- 
cussion), and  non -viscous.  Since  baffles  are  not  considered  and  the  tanks 
are  relatively  large,  viscosity  has  little  effect  on  the  fluid  motions. 

All  tank  motions , ( except  those  normal  to  the  mean  free  surface ) , are 
restricted  to  small  accelerations.  In  addition,  the  angular  motions  of  the 
tank  are  restricted  to  small  displacements.  The  free  surface  displacement 
and  slope  must  be  small.  These  conditions  insure  that  the  equations  can 
be  linearized.  This  is  consistent  with  the  rhirI 1 perturbation  methods  used 
in  airplane  dynamics  and  does  not  necessarily  limit  the  value  of  the  results, 
(in  this  connection  it  should  be  noted  that  the  above  restrictions  are  suf- 
ficient for  linearization,  but  not  all  necessary.  The  necessary  and  suffi- 
cient condition  is  probably  that  free  surface  slopes  and  displacements  be 
small  when  measured  relative  to  axes  fixed  in  the  tank. ) 

Only  rectangular  tanka  are  considered  here  in  detail.  The  cylin- 
drical tank  of  circular  cross-section  with  free  surface  normal  to  the  axis 
has  been  treated  in  Ref.  4. 

Free  Oso-j  nations 

To  determine  the  nature  of  the  free  oscillations  in  a stationary 
tank  the  procedure  is  as  follows.  Since  the  fluid  is  assumed  non -viscous 
and  incompressible,  the  velocity  potential  is  chosen  to  satisfy  the  Laplace 
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equation*  At  the  vails  the  normal  velocity  components  must  he  zero*  On 
the  free  surface  the  pressure  must  be  zero,  and  in  order  to  satisfy  this 
condition  vith  linear  functions  it  is  necessary  to  assume  mm>n  displace - 
■ents  of  the  free  surface.  For  rectangular  tanks  this  problem  is  readily 
solved. 

Forced  Oscillations 

In  studying  the  forced  oscillations  of  the  fluid  to  determine  an 
equivalent  mechanical  system  it  is  sufficient  to  study  the  response  of  the 
fluid  to  simple  harmonic  motions  of  all  frequencies*  If  the  system  re- 
sponds correctly  for  each  frequency  it  vill  respond  correctly  to  any  combi- 
nation of  these  frequencies,  since  the  problem  has  been  linearized.  An 
arbitrary  forcing  function  can  be  constructed  by  superposition,  so  the 
mechanical  system  is  equivalent  to  the  fluid  for  any  arbitrary  disturbance* 

The  forced  oscillation  requires  new  boundary  conditions  corresponding 
to  the  motion  of  the  tank  vails*  Comparatively  simple  velocity  potentials 
can  be  constructed  to  satisfy  these  conditions.  However,  these  potentials 
create  a pressure  variation  at  the  free  surface  vhlch  cannot  actually  exist. 
It  is  then  necessary  to  superimpose  additional  potentials  having  the  time 
frequency  of  the  forcing  function  and  the  space  character  of  the  free  oscil- 
lations to  cancel  out  this  unwanted,  pressure  variation.  The  resulting 
potential  satisfies  the  boundary  conditions  for  the  moving  tank  and  complies 
vith  the  zero  pressure  condition  on  the  free  surface. 

Tank  Motions 

Any  tank  motion  of  the  type  considered  here  can  be  constructed  from 
motions  of  translation  along  the  three  axes  in  space  and  a unique  set  of 


rotations  about  the  three  axes.  (One  axis  is  taken  as  perpendicular  to  the 
mean  free  6urf ace . ) Motions  normal  to  the  free  surface  vill  be  discussed 
later.  This  leaves  three  basic  types  of  motion  to  be  considered  at  present; 
translation  parallel  to  the  free  surface,  rotation  about  an  axis  parallel 
to  the  free  surface  and  rotation  about  an  axis  perpendicular  to  the  free 
surface.  For  convenience  these  motions  are  referred  to  as  horizontal, 
pitching  and  yaving.  The  horizontal  motion  is  comparatively  easy  to  analyze, 
but  the  pitching  and  yaving  cases  are  considerably  more  complex.  Ihe  yaving 
case  requires  the  use  of  a 3-dimensional  solution  of  Laplace's  equation  end 
introduces  a nev  set  of  natural  frequencies.  The  pitching  and  horizontal 
motion  cases  require  only  tvo -dimensional  solutions  and  possess  a common  set 
of  natural  frequencies . 

Alternative  Methods 

The  original  method  used  for  solving  the  forced  oscillation  problem 
made  it  possible  to  set  up  equivalent  mechanical  systems,  and  gave  expressions 
for  forces  and  moments  only  as  functions  of  the  forcing  frequency.  If  it  is 
not  desired  to  use  the  mechanical  analogy,  an  alternative  process  is  to  set 
up  the  entire  boundary  value  problem  in  terms  of  the  Laplace  transform  of 
the  velocity  potential  (instead  of  using  the  velocity  potential  itself). 

Bii6  leads  to  expressions  for  the  Laplace  transforms  of  the  forces  and  mo- 
ments, vhich  are  included  in  this  report.  If  Laplace  transform  methods  are 
used  for  solving  the  complete  dynamics  problem,  then  these  force  and  moment 
transforms  enter  directly  as  terms  in  the  equations. 
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DEVELOPMENT  OF  EQUATIONS 


Boundary  Conditions 

In  the  original  development  of  the  equations  for  the  moving  tank  it 
vas  assumed  that  the  boundary  condition  at  any  vail  would  be  satisfied  only 
at  the  mean  position  of  the  vail.  In  the  case  of  horizontal  tank  motions 
such  an  assumption  implies  severe  restrictions  on  the  motions.  Since  these 
restrictions  are  both  undesirable  and  unnecessary,  they  have  been  eliminated 
(for  horizontal  motion)  in  the  following  development. 


Consider  a rectangular  tank  filled  to  a depth  h with  fuel.  Let  the 
vails  and  bottom  of  the  tank  be  given,  in  irrotational  coordinates  (x,y,z) 
which  follow  the  horizontal  motion  of  the  tank,  by  x = - ^,  y = and 

* ■ - 7j.  If  the  horizontal  motion  of  the  tank  center  is  given  in  inertial 


coordinates  (x,y,z)  by 


x =XM 

y = Y( t) 


then  the  moving  coordinates  are  related  to  the  inertial  coordinates  by 


•X  = X -X(t) 


z - z 


Thus  a function  f(x,y,z,t)  in  the  inertial  coordinates  would  transform  in 
the  moving  coordinates  to 


(3) 
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The  velocity  potential  ft  (x,y,z,t)  which  describes  motion  of  the  fuel  in 
response  to  tank  motion  must  satisfy  the  conditions  that  at  any  fluid  boundary 
the  normal  velocity  of  the  fluid,  - ^ , be  equal  to  the  velocity  of  the 
boundary  normal  to  itself*  If  the  velocity  potential  is  expressed  in  the 
form 


<P(*fyz,t)  = -X(6)x-  + <P  (X,y,z,*) 


(M 


then  <P 1 must  satisfy  the  conditions  that  at  any  fluid  boundary  the  normal 
velocity  of  the  fluid  relative  to  the  moving  coordinates,  - ft)J^  , be  equal 
to  the  normal  velocity  of  the  boundary  relative  to  the  moving  coordinates* 

In  the  case  of  pure  horizontal  motion,  ft)  /=  O on  the  tank  walls  and  bottom. 
If  motions  in  pitch  and  yaw  involve  only  small  angular  displacements  of  the 
tank  from  its  original  orientation,  then  the  boundary  conditions  can  be 
applied  with  sufficient  accuracy  at  x » - p y = - ^ and  * « - 

If  all  tank  motions  are  restricted  so  that  they  involve  accelerations 
which  are  small  compared  to  the  acceleration  due  to  gravity,  and,  in  the 
case  of  sinusoidal  motion,  the  forcing  frequency  is  not  too  near  the 

frequency  of  a mode  of  free  oscillation  of  the  fuel,  then  ftg  , ft-  , and 

— f — 

<p~  will  be  small.  Also  the  slope  of  the  free  surface,  > wil1 

small. 

Neglecting  the  squares  and  cross  products  of  these  small  quantities, 
the  linearized  free  surface  conditions  in  the  inertial  coordinates  are^ 


(5) 


where  F(t)  1b  an  undetermined  function  of  time.  By  Eqs.  (2)  and  (3) 

so  that  the  free  surface  conditions  become,  in  the  moving  coordinates: 


$Pl  = ft  + Y*+r(t) 


(7) 


(x>y>7h>  (8) 

The  quantity  -L-yr^L  in  Eq.  (7)  can  be  set  equal  to  zero,  since 

F(t)  is  an  arbitrary  function  of  time,  and  and  Y are  functions  only  of 
time.  This  quantity  does  not  affect  the  fluid  motions  since  it  introduces 
no  pressure  gradients.  Eliminating  ^ between  Eq.  (7)  and  (8)  the  combined 
condition 

#tt  (x-V’T-*)  *7  4 (*->#>*) = ° <9) 

is  obtained. 

If  it  is  assumed  that  the  horizontally  oscillating  tank  is  equivalent 
to  a stationary  tank  with  oscillating  horizontal  force  field,  the  preceding 
conclusions  are  reached  intuitively. 

To  summarize,  the  problem  now  reduces  to  determining  <p  (x,y,z,t)  to 


satisfy  the  following  equations: 


f Ax  &t**) + fyNb**) + = 0 

(at  tank  boundary)  =»  normal  Telocity  of  tank  boundary 

At  (x>y>£’*)  + (x>y4,*) = ° 

In  the  following  derivations,  three  of  the  five  basic  tank  motions 
are  considered.  They  are:  horizontal  motion  parallel  to  the  x-axis, 
pitching  motion  about  the  y-axis,  and  yawing  motion  about  the  z-axls.  Solu- 
tions for  horizontal  motion  parallel  to  the  y-axis  and  pitching  motion  about 
the  x-axis  can  be  obtained  by  interchanging  a and  b in  the  first  two  types 
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of  motion  considered. 

Sinusoidal  motions  are  considered  first.  Arbitrary  tank  motions  are 
treated  by  use  of  Laplace  transforms.  In  each  case  it  is  assumed  that  the 
motion  involves  accelerations  which  are  small  in  comparison  to  the  accelera- 
tion due  to  gravity  and  that  the  frequency  of  sinusoidal  motion  does  not 
have  a value  in  the  Immediate  neighborhood  of  the  frequency  of  a natural 
mode  of  free  surface  oscillation  in  a stationary  tank.  Motions  in  pitch 
nrvi  yaw  are  assumed  to  involve  only  small  displacements  while  horizontal 
motions  may  involve  large  displacements  and  velocities.  This  is  consistent 
with  the  linearized  treatment  of  aircraft  dynamic  problems.  Horizontal  dis- 
placements of  several  tank  widths  may  occur  during  side  slip  or  in  the  case 
of  a tank  located  far  from  the  center  of  gravity  of  an  aircraft  undergoing 
small  yawing  motions. 

Sinusoid Horizontal  Motion  Parallel  to  the  x-Axis 

If  the  horizontal  motion  of  the  tank  is  given  by 


JC  (t)  = As&hzcvt 


(10) 


then  the  velocity  potential  <f>  (x,y,t)  must  satisfy 

- 0X  (i-r,z,-6)  =Acoa*scot  (n) 

<t>z  (x,-4’*)‘°  (12) 

The  velocity  potential  which  satisfies  Laplace ' s equation  and  Eqs.  (ll)  and 
(12)  is 


<j>  - -A  co  c&s  cut  jx 

+ ^/B7} s4m[(2n+i)% ‘xJcGsJifCz 71+/)tz (% +y)]\  ^13 * 

7Z=0  J 

The  term  proportional  to  x satisfies  the  velocity  condition  on  the  tank  vails, 
while  the  terms  in  the  summations  are  needed  to  satisfy  free  surface  conditions. 
Applying  the  free  surface  condition  to  Eq.  (13)  gives  the  identity 

Au >c&j  clt&Jx  -h  Zn$^urJi^ntf { &-] 

+T/Bn  fco&c<!#A[(2n+i)Vzt]  ~%(z  n+i)^^Ln/i jjfen+iflY  'af]') 

n'-°  r „,i  <nt) 

~){.4in[&>»-0'3rJ >=o 

Expanding  x in  a Fourier  series  and  then  equating  coefficients  to  zero  in 
the  Identity,  one  obtains 


B i (-ft  4 <*■  — _ 


j 
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By^  becomes  infinite  if 


U)2=  O .=  p (Z  n -h) ~-bmh[(2n+i)  IT ^ ] 


CO^  is  the  frequency  of  the  nth  odd  harmonic  of  the  fundamental  mode  of 
free  surface  oscillation  in  a stationary  tank.  The  velocity  potential  is 


4 = 


v/  ^<iun  [(Zn+ijH x]c<xdi  ((2n+i)zu(2  +iBj  v 
o&zA[(2  n+i)7T-bc.J  J 

The  elevation  of  the  free  surface  above  z = - obtained  with  the  aid  of 
E(J*  (7)  OO  , 

? = A cot 5 jx  { 'a&s) 

V (19) 


The  horizontal  force  exerted  by  the  fluid  on  the  tank  side  vails  is  obtained 


by  integrating  the  difference  in  pressure  at  equal  elevations  on  vails 


- AuAdiwcvt 


\/  (a?b  SficSb  fb  2a.tm/iknH)1nll . 9-  If  CO2  \) 

Xfe  is  <znX)r  J+4Jwsrjr 


Sinusoidal  Pitching  About  The  y-Axis 

If  the  pitching  motion  of  the  tank  Is  given  by 


Q ~ JB  CO  ~t 


where  B is  the  maximum  angular  displacement  of  a tank  vail  from  its  mean 
position,  then  the  velocity  potential  must  satisfy  the  boundary  conditions: 


=+JBcuz  c&Scu-£ 
~<pz  ~ ~Bcux  c&dcot 


The  velocity  potential  which  satisfies  these  boundary  conditions  is 


r co 

<f>  = - B ooc&ioot ' 

171=0 


V—+ 

7r3(zm-if 


sdimA,  jfen+f) 


+ £a»  ax  P^x]cmA[ &Z(z +£)] 


(27) 


Applying  the  free  surface  condition,  Eq.  (9)/  the  identity 


r 2r^  4 k 2 [(£  n+fijC x] 
| ^Zj  Tf^n-HfaxUi [fcmDffi-] 


+ c gfi_  a77  ^ a?^n  [&*+')  ZU*] 

^ J(2  n-tf)  3aBrA  [(27j+/}§;hj 


71=0 


+ 


YjA^  (ofa*##*]  inA[zn£h]') 

yCcozfenrrS;] 
lL/Bn  (^c&dhji&i+i) JF ^~^(27Hi)^2uAA^2n+i)S^f) 


(28) 


is  obtained. 

Expansion  of  the  first  summation  in  Eq.  (28)  in  a Fourier  series  gives 

^Ix^oriJi  fcnHtfx] 


vp  ^Ltx_^22BJzL——L K—d 

£o  VSM3*^fe*+')g£] 


oo 

-2 


(-1)4, 


Z/  7rz(2n+i)z  [£  (2n.H)V' 

7 1=0  u 


(29) 


IAAnlfcn+/)%ij 


which  upon  substituting  into  Eq.  (28)  and  equating  coefficients  to  zero 
gives  after  simplification 


A 


/ 

71 


= o 


(30) 


fa  AcunA 

(2n+j)W 
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The  velocity  potential  Is 


(fi= - JBcuc&jcot* 


--  4 

C } ~3'en+ff[  C04Aj(zn+/)jn±-] 


rr3(sn+i) 

a*  Jin  [&w)zx]c^[(2n+Q%Cz--T)]  1 

MwA j(27i+/)]fh.]  3 


(32) 


YjH) 

71=0  V 


CO 


fk  2a1amA[(2w)!£] 


Jr  ** 


+ 


tk] 


x 


CoaAjfcn+OSrA] 


The  elevation  of  the  free  surface  la  obtained  by  use  of  Eqs.  (7)  and  (29)  • 


7*hp.  horizontal  force  exerted  by  the  fluid  on  the  tank  vallB  la  obtained  by 
use  of  Eq,  (20) 


The  total  moment  about  the  center  of  the  tank  la  obtained  from  Eq,  (22)  where  In 
the  second  integral  (g  h)is  replaced  by  g ( fl-BX^neoi:)  to  take  Into  account 
the  unequal  depth  due  to  rotation  of  the  tank  through  the  angle  B 4^71  CO  i; . 
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[ ^77^(2 n+ 1)4 (2n+/)7T  co^ 


-L 


f'  8P  h3b  fa,  3/1  'Mi[(Znti)^x]  7 a &p*3b 

^7r+(2n+Q+Lz  (Zk+Ott  J /-s$- 7r*(zn+/)4 


k,  2 cl  tcorih[(Z  ti  +i)z  J f- 
Z"  (2n+i)ir  ; 


■^5*  -4^n  tu  -b  ^jjr 


Sinusoidal  Yawing  About  The  z-Axls 

If  the  yawing  motion  about  the  z-axis  is  given  by 


y - C-^bn  CaJ  -t  ^ 

where  C is  the  maximum  angular  displacement  of  the  tank  walls  from  their 
mean  position,  then  the  velocity  potential  must  satisfy  the  boundary 


conditions: 


-0x(~jr>y->  z>t)  ~ -Ccuy.  c&d  co  -t 

~<Py(x>~4~’  z>*)  = ^ t 


(39) 


r 


MT  = BcJAun. wt!  f^^dA  fh  a.tvnk[an+i)ik]  _± 

) ^Tr+(2mifl7  (zn+i)7T  col \j 


V'  8P  h3b  [a,  3k  Mfcn+Ott]  h 

£/oTr+(zn+0*  Lz  (Zn+Oir  J 7r*(zn+0 


(Zm-H)7T  J ir*(Zh+0+ 


yf/l  2 a.  tcto%k[(zn+f)^k:]  %-  If  OJ2 

A[*  (2n+t)rr  **?/  [co£-cc 
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Slnusoldal  Yawing  About  The  z-Axls 

If  the  yawing  motion  about  the  z-axla  is  given  by 


'P  = C^dPri  OJ  ~t  (3^ 

where  C is  the  maximum  angular  displacement  of  the  tank  walls  from  their 
mean  position,  then  the  velocity  potential  muBt  satisfy  the  boundary 


conditions: 


— (J)^ pt. ~ -'Qcu'ij,  c&d  u>-6 
-<fL(x,±-§-,z,t)  = Cujx  cva  cut 


4 ° 


and  the  Laplace  equation  In  three  dimensions.  The  velocity  potential 
satisfying  these  conditions  is 


j r „ , f fV  A*  * 

<(>  - I 

_ a? Ain  [(2n+j)  1 

&*A[(Z  n h)  g£]  J 

J (4o 

oo  oo 

Aim[(z  7i+i)gx]Amfe7n+f)jy] 

771^0  n-o  » L 

X c<^[^V6Y^f+^^‘f(z  +£jj | 


where  terms  containing  even  multiples  of  n or  m have  been  omitted  (since 
they  contribute  nothing  to  the  velocity  potential) . The  term  involving 
only  the  single  elimination  is  the  part  of  the  velocity  potential  which  satisfies 
the  velocity  conditions  on  the  tank  boundaries,  while  the  double  summation 
term  is  needed  to  satisfy  free  surface  conditions. 

Applying  the  free  surface  condition  to  Eq.  (kO),  the  following 
identity  is  obtained: 


~Cod  C&4  cO  t 


sfS  ,i”v  4-b2Ain/i  l&m+nX]  ■ n 


">frJ 


2 ,.n  4a 

" U ’ TT^nHfc^^)^] 


(Equation  continued  on  next  page) 


Upon  substitution  of  the  Fourier  expansions 


s<Urdi [feTn  +dTX]  - 4ab'\' (~pA^n [(zn+l)  2X] 
(2™+i)c<x^lfe-rn+/)jfr]  ^ n=o  [6£(2n+/)i+a2(2m+i)z] 

^7zy/;^^y/;^7  ^rn^o 


into  Eq.  (lH),  the  coefficients  A™  n are  found  after  simplification  to 

J v 


^7^,7?  -An ,n  [co^n-coy 

£ (~iyi^/6ab[b%ni-i)z-<^'m-i-<) ]mcA [^^bZ(2n+ bZ+^2(2m+,f^ 

7T*(2m +tf(z n+t)z  [b 2 (2 n h)2+  a2(2m+if] 

where 

C0m>n= fStib2(zn*f+a‘(im.^^‘[&  W&n+O**^™^ h]  <45) 
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The  elevation  of  the  free  surface  is  found  with  the  aid  of  Eqs.  (7),  (42) 
and  (43). 


- r af  J V ^(-i)m+nJ6a6l^K+dZ-afe™  nf] 

_co6co  | Zj2-/  Tri{em-hif(2n  +i)%2(2n -tif+dfem+i)2] 


Since  the  notion  iB  purely  yawing,  the  fuel  exerts  no  net  horizontal  force 
or  manent  about  the  x or  y-axis.  The  moment  about  the  z-axis  or  yawing 
moment  exerted  by  the  fuel  on  the  tank  is  obtained  frcm  the  following 
expression 


Jx 

'2 


A 

'Z 


Mz  = -\P 


A 

2 


A 

2 


a, 

'2 


a 

2 


[fy  (x}  2 J 2lZ)t)JxdX>dz 


(47) 


Substituting  ^rnn  int0  EtJ*  (4°)  and  aPP1yin€  EQ*  (47)  the  yawing  moment 
is  found  to  be  r 

+b4UnA[&w)ffr]J-7ip*  bh  (az+ &) 

y-v  b3 [b^tin+zf- a? fem+if]  f coz  \ 

m-o  7t -o  %■  It 8 [b Z (2 n +tfi“  ^C2- 171  (zn+i)  (pi+0  i^7.n  °^J 


(48) 
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The  first  summation  of  the  right  side  of  Eq*  (48)  is  the  effective  moment 
of  inertia  of  the  fuel  about  the  z-axis. 

Arbitrary  Horizontal  Tank  Motion 

The  Laplace  transform  (vith  respect  to  time)  of  the  boundary  conditions 
for  arbitrary  horizontal  tank  motion  is 


(t£,z,i)}  = - 


£{<t>z(x,-b^}-°  (50) 

where  ie  the  displacement  of  the  tank  in  inertial  coordinates.  Under 
the  assumption  that  up  until  the  time  t = o the  free  surface  is  undisturbed, 
the  Laplace  transform  of  the  free  surface  condition  is 


jb'Xft  &,£,*)}  = 0 (51) 

The  Laplace  transform  which  satisfies  the  conditions  given  by  Eq.  (it-9) ^ (50) 
and  (51)  is 


X{4}  = -X{X(t)} 

V [ ^ 4. CL  / i>Z  \ 


The  Laplace  transform  of  the  elevation  of  the  free  surface  above  z = ^ 1b 


given  by 


The  Laplace  transform  of  the  horizontal  force  exerted  hy  the  fluid  on  the 
tank  walls  is  given  by 

Xfa*)}  = *[£{$&>  z‘*)}-£{$(~f’  E,-i)}]dZ  (55) 

°~T 

xffi*-  ?£{*■&} f°aM  -£  (56) 


The  Laplace  transform  of  the  total  moment  exerted  by  the  fluid  on  the  tank  is 

A 

X{Mrfc)J  =pbf>  j % Zjtylzd z 

U~7T  Pf  , (57) 

+pbf> bCfyfc  x dx 


&{MT}=-j>£{X(t)} 


Y fp  C&b  ST'CCfc.  0/Oajb/h  2atamhfen+i)£-thl  §\f  £f_ 

A//  -/IT  L>0  f Tr*f>nHf\2  czn+i)  n 


The  Laplace  transforms  of  velocity  potential,  free  surface  elevation, 
horizontal  force  and  total  moment  representing  the  response  of  the  fuel 
to  arbitrary  pitching  and  yawing  motions  of  small  amplitude  are  given  below 
without  proof. 


DOUGi 


Arbitrary  Pitching  Motion 


/ji/mfa-n+OjiZ 'J'&nJifen+Ox-xj 


fi ) y>  (a  m)  rO (tO  ± A j^t"ULL^n 

^ i tj  _c^(  JLn^o  c^fij(2n-hi)-^^j 

, a.*AML F(2n-H)7rX Icozklfen+ijfcfZ - jtfj)  ST',_f}*  ±a,  / AZ 

L T,  nr.  TArTn  / 2u  irfy-n+i^X^+crf. 


4UiA[(zn+i)£kJ 


Aj0  irfyn+if^+co 


'h  20  ia^A[^n-hl)j- j] , %-  \^An/^nH)^C<^[(Zn^)a(Z^^ 


y f k 20  iaoihLl 

^ ( 2 (2n+i 


c&nJiffzn+^/ij 


x/Mfe 


£{$}•'&&&} 

\/ FncJb  V'c0n  epa?b  fh  2a.&nA[(inH)zitl  $-)(  ’P 

*/'  iz  aU  ? * (Zn-ht)Tr  a%)\p-H. 


r+v, 


S^r  ) / 'i=>fA/j-\'\  r ^ 8pa?b  fh  a J+-2- 

(ZnH)TT 


A~\  ffph3 b fa,  3 A tamJiffen+i)  z A,  J 


&pALA  — 

-jf^o  T4(Zn-n)4\Z 


fern-Off 


(Equation  continued  on  next  page) 


Y(t)\ fc  a bffizn+i)*-  a*(zm+if] 

l i f L J \rn=o  n:o 


•where 


_-of 'cAn  e+p°?b*[t^f-^>n+if]  J_£l_)j 
nwnlo  & V9 [6fynH)*+tf2m-Hffon+i  faH)  \ +U^*>1j\ 


& >*1,71  is  given  by  Eq.  (44)  and  yJCi:)  is  the  angular  displacement 


of  the  tank  walls  in  yaw. 


EQUIVALENT  MECHANICAL  SYSTEM 


Consider  a rectangular  tank  of  length  a,  vidth  b,  and  mean  fuel 
depth  h undergoing  small  amplitude  sinusoidal  oscillations  of  frequency  GO 
about  a set  of  fixed  axes.  Let  the  fixed  axes  coincide  vith  the  principal 
axes  of  the  fuel  when  the  tank  is  at  reBt  in  an  upright  position.  The 
positive  directions  of  these  axes  are  shown  in  Fig.  1.  Positive  forces, 
moments,  displacements,  and  rotations  are  defined  according  to  the  con- 
ventions of  a right  hand  system. 

The  problem  of  replacing  the  fuel  by  an  equivalent  mechanical  system 
resolves  itself  into  finding  a mechanical  system  for  each  of  three  types 
of  motion: 

(A)  horizontal  motion  parallel  to  the  x-axis  and  pitching  about  the 

y-axi8. 

(B)  horizontal  motion  parallel  to  the  y-axis  and  pitching  about  the 

x-axis. 

(C)  yawing  motion  about  the  z-axis. 

In  order  for  a mechanical  Bystem  to  be  equivalent  to  the  fuel  with  respect 
to  one  of  the  above  types  of  tank  motion,  it  must  exert  the  same  force  and 
total  moment  on  the  tank  as  the  fuel  does  for  that  type  of  tank  motion. 

Motions  of  type  (A)  and  (B)  are  of  the  same  character  so  that  it  is  only 
necessary  to  discuss  one  of  them. 

Motions  of  Type  (A) 

Consider  a mechanical  system  composed  of  a fixed  mass  M and  an  infinite 
set  of  undamped  Bpring-masses  j"  80  constrained  as  to  move  only  parallel 

to  the  bottom  of  the  tank  and  the  xz-plane.  Let  the  nth  spring -mass  have  a 
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spring  constant  suc^  that 

** » 

where  i8  frequency  of  an  odd  harmonic  of  fundamental  free  fuel 
oscillation.  Locate  the  nth  spring  masB  at  a vertical  distance  zn  and 
the  fixed  mass  at  a vertical  distance  Z as  shown  in  Fig.  2.  The  moving 
masses  may  be  thought  of  as  point  masses  while  the  fixed  mass  has  a distri- 
bution such  that  its  moment  of  inertia  about  an  axis  parallel  to  the  y-axis 
and  passing  through  2 « Z is  . 

If  the  tank  undergoes  a horizontal  oscillation  of  amplitude  A and 
frequency  00  parallel  to  the  x-axis,  then  the  equations  of  motion  of  the 
moving  masses  are  given  by 

= - k-n  (x^AMnoot)  (66) 

vhere  is  the  displacement  of  the  moving  masses  parallel  to  the  x-axis  • 

If  the  tank  undergoes  a pitching  oscillation  of  angular  amplitude  B and 
frequency  CO  9 the  equations  of  motion  are 

TnnXn=-kn(xn-S2„Mncot) ydn  cot  (67) 

The  horizontal  force  and  total  moment  exerted  by  the  mechanical 
system  were  obtained  with  the  aid  of  the  solutions  of  Eqs.  (66)  and  (67). 

A comparison  of  these  forces  and  moments  (designated  by  primes)  with  the 
corresponding  forces  and  moments  exerted  by  the  fuel  on  the  tank  is  shown 


below.  The  forces  and  moments  have  been  reduced  to  non-dimensional  form 


in  which 


Mp  = f>  abh  =>  total  fuel  mass  • 

**  g P abh  n total  fuel  weight 
g » vertical  acceleration 

I?  » effective  moment  of  inertia  about  the  y-axis  of  the  fuel 

y 

r^  « — * tank  aspect  ratio 

= = (.Zni-i)  TrA,tamJi[(Zn+i)irA,J 

rZ  kcos 


Sinusoidal  horizontal  motion: 


Fh__A 

WF~  h.' 


’McmAiSma,t \=L, 

Wfh  h J \/2/li  nA  Tf3(zn->-i)3/il 


v /)  fanhl&n + 0 . / ) l 

A p otn+ofstr 


Sinusoidal  pitching  motion: 


oo 

vPMa  i 


ffi- 


UOUJLUUX  yj  X wxv^ii 

V f / \ J_ 

(znH)f/l;  fnjjy^2 


If  the  forces  and  moments  exerted  by  the  mechanical  system  are  to  be 
Identically  equal  to  the  forces  and  moments  exerted  by  the  fuel,  then  the 
characteristics  of  the  mechanical  system  must  be  given  by 


7772k _ 8tariA[(2ni-i)rr/iJ 
Mf  ~ 7f3 (2n+i)3 A/ r 


z n _ / tcurvh^n-tiyj^i] 

k 2 (zn+i)3£fli 


(77) 
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Tty  _ Xfqj,  M_  (T\Z  FFlix  fzn\Z 
Mrh*  ~MFhz  Mf  [7TJ  £/0Mf\  fl  j 


(So) 


kkyv  = 8 tarJi  ffeni-Qlfc,]  (8l) 

Wf  ft*  (2n+i')z 

It  1b  Interesting  to  note  that  the  rigid  mechanical  system  has  the 
same  moment  of  inertia  about  the  y-axis  as  the  fuel  if  the  free  surface  vere 
constrained  by  a tank  boundary.  This  effective  moment  of  inertia  of  the 
fuel  is  found  to  be 


vhere  is  the  moment  of  inertia  about  the  y-axis  if  the  fuel  vere 

solidified. 


Approximate  expressions, 
Ifa^I  I and  Z/A 


good  to  five  significant  figures,  for 
which  lend  themselves  to  easy  calculation 


are  given  below. 
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(83) 


M 

Mc 


- /_  Qi35J?0/z  [tanA  D'/I,  + o. 05/80  J 


(84) 


for  Ay  - £ 


M 


= /.032049A,\p2yrJljj^i  + 0.05/80oJ 


m 


for 


j—  + (—  - OJM2A6  /"  / , / L, 

f^)  «<(^  h*-'**"  aWaffi*,  7 


z.^  ^ 

X 

The  mass  ratios  9ri0^MF  , 7/2^//] F , and  A/1  f/15 F have  been  computed 
from  Eqs.  (76),  (84)  and  (85)  and  are  plotted  as  a function  of  ri  in  Fig.  3. 
The  moving  masses  associated  vith  the  higher  harmonics  are  too  n«m i to  plot 
except  very  near  r^  = 0. . The  ratio  of  to  77L0  is  approximately 

2 o 

l/ (2n+l)  for  small  r^  and  l/(2n+l)  for  large  r-^%  These  ratios  axe  plotted 
on  a logarithmic  scale  verses  the  ratioy^/^  in  Fig.  4. 

The  ratios  of  the  arms  5q  ^ 2/  ^ and  Z to  have  been  computed  from 
Eqs.  (77)  and  (86)  and  are  plotted  in  Fig.  5.  As  the  tank  becomes  shallow, 
the  ratio  becomes  infinite,  but  Z.  approaches  0.46o6l3a  as  can  be  seen 

by  multiplying  Eq.  (86)  by  ar^  and.  letting  r^_  approach  zero. 
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The  moment  of  inertia  ratios  IF^/ 1$^  and  ly/ 1$^  have  been 
calculated  from  Eqs.  (80)  and  (83)  and  are  plotted  in  Fig.  6.  It  should 
be  noted  that  values  of  Ip  j -7^.  for  r^>  1 can  be  obtained  by  replacing 


W and  IJh 


rj_  by  l/r^  in  Eq.  (83). 


r y 


Fig.  7 is  a plot  of  the  non-dimensional  form  of  the  spring  constants 


, and  as  calculated  from  Eq.  (8l). 


In  order  to  get  a dimensionless  parameter,  the  horizontal  force  pro- 
duced by  the  mechanical  system  is  divided  by  the  force  which  would  be  produced 
if  the  fuel  were  solidified.  Fig.  8 is  a plot  of  this  ratio  as  a function  of 
frequency  for  a tank  aspect  ratio  of  0.25 . At  resonant  frequencies  the  force 
ratio  would  not  actually  become  infinite.  In  Ref.  3c,  Lamb  discusses  a limit- 
ing value  of  the  ratio  of  wave  amplitude  to  wave  length  which  permits  compliance 
with  free  surface  boundary  conditions,  nils  limiting  ratio  was  found  to  be 
approximately  0.1b2  beyond  which  one  may  assume  that  energy  dissipation  through 
splashing  would  occur.  Corresponding  to  any  amplitude  limitation  of  any 
particular  mode  of  fuel  oscillation  there  is  a force  limitation  in  that  mode. 

At  resonant  frequencies  for  the  tank  harmonics  the  maximum  force 
that  can  be  produced  in  the  nth  odd  harmonic  mode  is  l/(2n+l)2  times  the 
may i mum  force  that  can  be  produced  in  the  fundamental  mode  (n=o).  This 
gives  sane  Justification  for  discarding  the  moving  masses  associated  with 
modes  of  fuel  oscillation  beyond  the  fundamental. 
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Motions  of  Type  (C) 

Consider  a mechanical  system  composed  of  a fixed  mass  vith  mcment  of 
inertia  IZ  about  the  z-axis  and  an  infinite  set  of  moving  masses  constrained 
to  pivot  about  the  z-axis.  If  the  moving  masses  have  moments  of  inertia 
and  are  attached  to  torsional  springs  vith  spring  constants  ,n,  then 
the  equations  of  motion  of  the  moving  masses  are 


I rpr  - -sk  -c Ain  cot) 

r?ri,n  on  n / 


(87) 


771,7?  ^ 771,71 


vhere  Ca&TI  cu  ~t  is  the  angular  displacement  of  the  tank  in  yav.  If  the 
spring  constants  are  related  to  the  natural  frequencies  of  free  fuel  oscil- 
lation by 


^?7?, 7i  ' ^?77; n 


(88) 


then  the  yaving  moment  exerted  by  the  mechanical  system  on  the  tank  can  be 
found  from  Eq..  (87)  and  compared  to  the  yaving  moment  exerted  by  the  fuel. 
This  comparison  is  made  be  lav  in  vhich  the  primed  moment  is  the  moment  pro- 
duced by  the  mechanical  system  and  vhere 

I*_  m moment  of  inertia  about  the  z-axis  of  the  solidified  fuel 
T ■ effective  mcment  of  inertia  of  the  fuel 

Z J, 

/r  . .xk  = tank  aspect  ratio  in  the  xz -plane 
/0!  & 


/t,2 
i 

'3 


« tank  aspect  ratio  in  the  yz -plane 
P 

Az  Q 

a7  T 


fl , « Tta  d 0 tank  aspect  ratio  in  the  xy-plane 


^ ^7rv,7b  _ (Zm  7t  +1)*  fy/rih  (zm+i)  +(zn.+l)  J 
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SC  sz 


C a,  U)  A/^n  nT)t\ 

f 


J-%  ikL1**  Jti&T*  ^)-/ 


(90) 


If  the  tvo  yawing  mounts  are  to  be  Identical  for  all  forcing 
frequencies,  then 

Im  n 768 ftjfemti)' *-fcn+ifJ*  . 

XH  " vjP  71h,J\/Jt£(ZT7t+i)!+(z  n+if^ 


Iz  _ If*  _ y y Iyn,n 

?SZ  ^z  'mZ'oTvo  ZsB 


(92) 


It  is  interesting  to  note  that  if  r^  “ (2n+l)/(2m+l),  then  ?7n,n  =0- 
In  the  case  of  a tank  with  square  planfono  (r^  ° 1),  the  fuel  associated  with 
the  frequency  ^ofO  oscillates  symmetrically  about  the  planes  x = -y  so  that 
no  net  force  or  moment  can  be  produced.  As  far  as  the  fuel  associated  with 
the  natural  frequency  ^ is  concerned,  the  planes  x » ±(2q+l)  a/2(2n+l), 

and  y = ±(2q+l)  b/2(2m+l)  where  q is  an  integer  could  be  taken  as  tank 
boundaries.  Thus  the  tank  could  be  divided  into  (2m+l)  times  (2n+l)  smaller 
sub -tanks.  If  r^  = (2n+l)/(2m+l),  the  smaller  sub-tanks  are  square  so  that 
no  net  forces  or  moments  are  produced  and  _Z"^  j^O  . 
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It  is  probably  sufficiently  accurate  to  consider  only  the  moment 
produced  by  the  moving  masses  with  moments  of  inertia  1 0,0 ^ Z/1ot  1 2,0, 

and  I0iz  . The  ratios  of  these  moment  of  inertia  to  ^ bave  been 
confuted  from  Eq.  (91)  as  a function  of  r^  for  the  shallow  tank  case,  r-|_=r2=0. 
They  are  shown  in  Figs.  9 and  10.  To  obtain  the  value  of  Tyri  yy/ for 
rj_  ^ °>  one  simply  multiplies  Zyyj  for  r^  = 0 by 

tank  [flfo/  '\jhz3(Zrn+i)z+  (zn+^] j 7 

To  obtain  values  for  r^  s'  1,  replace  r^  by  and  change  A,t  to  A. 2 . 

3 

The  ratio  Fz was  calculated  from  Eq.  (83). 

Mechanical  Systems  for  Combinations  of  Tank  Motions 

A single  mechanical  system,  that  will  respond  the  same  as  the  fluid 
to  all  tank  motions,  cannot  in  general  be  constructed  from  the  simple 
mechanical  systems  discussed  without  introducing  negative  moments  of  inertia . 
However,  it  is  always  possible  to  construct  a single  mechanical  system  that 
will  represent  the  fluid  response  to  both  yawing  motion  and  motion  of  type  (A) 
or  yawing  motion  and  motion  of  type  (B).  For,  consider  the  mechanical  system 
that  represents  motions  of  type  (A).  The  fixed  mass  can  be  distributed  in  the 
yz -plane  so  that  the  center  of  gravity  remains  fixed  and  its  moment  of  inertia 
is  Xz  about  the  z-axis  and  -Z"^.  about  an  axis  passing  through  B - Z.  and 
parallel  to  the  y-axis.  Tbe  moment  of  inertia  can  be  separated  into 

the  components  Zg  , 1 0}0,Z o,  / etc.  of  the  mechanical  system  representing 
the  response  of  the  fluid  to  yawing  motion. 
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DISCUSSION 

Seme  problems  arise  concerning  the  proper  use  of  the  preceding 
material,  and  a fev  of  these  are  discussed  below.  Also  Included  are  some 
topics  vhich  contribute  only  to  the  general  description  of  fluid  oscillations. 

Mechanical  Systems 

For  aircraft  haring  a vertical  plane  of  symmetry  through  the  longi- 
tudinal axis,  the  longitudinal  stability  and  the  lateral -directional  stability 
can  be  treated  independently. 

In  order  to  handle  the  longitudinal  stability  it  is  necessary  to  knew 
the  forces  and  moments  produced  by  translation  of  the  tank  and  by  pitching 
motions.  A single  mechanical  system  has  been  set  up  to  represent  the  tank 
for  both  of  these  motions. 

The  lateral -directional  case  is  more  complicated  since  in  general 
the  effects  of  rolling  motion,  translation  and  yawing  must  be  considered. 

A single  mechanical  system  has  been  devised  to  represent  the  fuel  tank  for 
all  of  these  motions,  provided  the  translation  is  confined  to  the  plane  of 
the  rolling  motion.  The  more  general  case  involving  also  longitudinal 
translation,  (vhich  would  appear  with  wing  tip  tanks),  cannot  in  general 
be  represented  by  a single  system  of  the  present  type,  unless  imaginary 
dimensions  are  introduced.  In  such  a case  different  systems  could  be  used 
to  obtain  forces  and  moments  in  different  planes. 

It  is  also  possible  to  use  the  Laplace  transforms  of  the  forces  and 
moments  directly,  without  regard  to  mechanical  systems,  if  the  transform 
method  is  being  used  for  the  complete  problem. 


In  some  cases,  particularly  where  tanks  are  located  far  from  the  C.G. 
of  the  airplane,  the  effect  of  angular  notions  of  the  tank  nay  be  negligible 
by  ccaqparison  with  the  effect  of  translation  of  the  tank.  In  such  cases 
sene  simplifications  nay  be  possible. 

Approximate  Mechanical  Systems 

The  forces  and  nonents  produced  by  the  higher  harmonics  of  the  fluid 
.notion  will  be  negligible  in  many  cases.  This  means  that  the  equivalent 
mechanical  system  can  frequently  be  assumed  to  consist  of  one  fixed  mass 
and  one  moving  bubs. 

If  the  airplane  has  a natural  frequency  quite  close  to  a higher  harmonic 
of  the  fluid  notion,  it  night  be  assuaed  that  it  would  be  necessary  to  retain 
the  moving  mass  corresponding  to  this  harmonic.  However,  consideration  of  the 
practical  limitation*  on  wave  heights  may  alter  the  picture.  These  limi- 
tations may  arise  from  the  effects  of  viscosity,  cavitation,  or  through  the 
presence  of  non-linearities  in  free  surface  boundary  conditions.  For  example, 
in  Ref.  3c  a limitation  on  the  ratio  of  wave  height  to  wave  length  is  suggested. 
Such  a "limit”  might  mark  the  Introduction  of  splashing  with  consequent  energy 
dissipation. 

Tanks  of  Different  Shape 

Tanks  of  other  than  rectangular  shape  could,  be  studied.  Stewart  and 

IL 

Lore 11  have  studied  the  forced  oscillations  in  & cylindrical  tank  of  circular 
cross-section  with  free  surface  normal  to  the  hxis.  Some  information  on  the 
free  oscillations  in  a few  other  tanks  of  simple  shape  is  given  in  the 
appendix  to  Ref.  6.  Very  few  cases  involving  variable  tank  depth  have  been 
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solved  exactly.  Although  there  are  mathematical  difficulties  in  getting 
exact  solutions  to  many  cases  which  appear  physically  simple,  it  is  still 
possible  to  get  approximate  theoretical  solutions  in  such  cases.  Also, 
it  should  be  comparatively  easy  to  determine  experimentally  the  fundamental 
frequency  for  a tank  of  arbitrary  shape. 

Saffle  Locations 

For  the  fundamental  mode  in  the  rectangular  tank  the  highest  vertical 
velocities  of  the  fluid  occur  at  the  ends  of  the  tank  in  the  free  surface. 
The  highest  horizontal  velocities  occur  in  the  free  surface  at  the  center 
of  the  tank.  These  results  are  easily  verified  from  the  expression  for 
velocity  potential. 

It  seems  probable  that  the  above  locations  are  best  for  baffles  since 
high  local  velocities  are  required  for  a high  rate  of  energy  dissipation. 

The  practical  problem  is  not  as  simple  as  this  would  suggest  however,  since 
it  may  be  necessary  for  the  baffles  to  be  effective  over  some  range  of  fuel 
heights  in  the  tank. 

Effective  Value  of  g 

If  the  fuel  tank  experiences  a large  constant  acceleration  normal  to 
the  mean  free  surface,  then  the  value  of  g is  effectively  changed.  If  the 
effective  value  of  g is  quadrupled,  then  the  natural  frequencies  of  fluid 
oscillation  are  doubled.  If  the  effective  value  of  g approaches  zero,  as 
in  a free  fall,  then  the  natural  frequencies  approach  zero.  Such  effects 
would  probably  be  of  interest  in  connection  with  missile  dynamics,  and  might 


be  of  importance  in  aircraft  maneuvers. 
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To  get  the  correct  equivalent  mechanical  systems  and  force  and  moment 
expressions  it  is  only  necessary  to  replace  g by  the  effective  value  of  g in 
the  formulas.  (This  results  in  the  introduction  of  an  effective  fuel  weight 
also.) 

Variable  Normal  Accelerations 

It  has  been  mentioned  that  a constant  acceleration  normal  to  the  mean 
free  surface  changes  the  effective  value  of  g and  alters  the  natural  frequencies 
of  the  fluid  in  the  tank.  It  can  also  be  shown  that  variable  normal  accelera- 
tions may  alter  the  amplitude  of  the  wave  motion.  In  particular  a periodic 
normal  acceleration  of  the  correct  frequency  may  feed  energy  into  an  existing 
wave  motion.  This  occurs  through  a non-linear  coupling,  and  the  equivalent 
mechanical  system  does  not  correctly  represent  the  fluid  to  this  order.  Such 
terms  are  not  ordinarily  considered  in  airplane  dynamics. 

Non-Linearities 

When  the  amplitude  of  the  fuel  motion  is  large  the  effect  of  non- 
linearity may  become  important.  The  non-linearity  is  associated  with  boundary 
conditions  at  the  free  surface,  not  with  the  Laplace  equation,  which  is  linear. 
The  mathematical  difficulties  introduced  by  the  non-linear  effects  are  too 
great  to  permit  investigation  here.  The  region  of  applicability  of  the 
linearized  equations  can  beBt  be  determined  by  ejq>erimento 

Viscosity 

The  effects  of  viscosity  are  small  for  unbaffled  tanks  of  the  size 
considered  in  aircraft.  If  baffles  were  present  this  would  certainly  not 
be  the  case,  since  the  actual  flow  pattern  can  be  radically  different  from 
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the  potential  flow  pattern  in  such  cases,  and  heavy  damping  may  he  introduced. 
For  this  reason  it  is  of  little  interest  to  investigate  the  potential  flow 
pattern  where  baffles  are  involved.  The  scope  of  the  report  is  therefore 
limited  primarily  to  determining  whether  or  not  baffles  may  be.  necessary. 

It  is  then  possible  to  suggest  locations  where  baffles  might  be  found  most 
effective.  This  is  based  entirely  on  the  velocity  distribution  in  the  un- 
baffled tank,  it  being  assumed  that  the  greatest  energy  dissipation  would 
be  produced  by  locating  baffles  in  the  regions  of  highest  velocity. 

Free  Surface  Shapes 

It  is  of  some  interest  that,  in  the  free  oscillations  of  fluid  in  a 
tank,  an  arbitrary  shape  of  the  free  surface  does  not  in  general  recur 
periodically.  For  example,  in  a stationary  rectangular  tank  if  the  free 
surface  is  initially  an  inclined  plane,  it  will  never  assume  this  shape 
again.  This  is  for  a non-viscous  fluid,  which  would  continue  its  oscillations 
indefinitely  without  loss  of  energy.  A sinusoidal  free  surface  whose  wave 
length  is  2/n  times  the  tank  length  (where  n is  an  integer)  will  repeat 
periodically.  However,  each  different  wave  length  repeats  in  a different 
period  of  time  and  these  periods  are  not  in  general  related  to  each  other 
through  rational  numbers . This  of  course  explains  the  non -repetition  of 
the  plane  free  surface  since  it  is  initially  composed  of  these  sinusoidal 
elements  (by  Fourier  analysis),  end  after  the  first  instant  they  never 
have  simultaneously  the  correct  relative  magnitudes  required  to  represent 
a plane. 

Superposition  Processes 

Knowledge  of  the  steady  state  response  of  the  fuel  to  all  frequencies 


of  oscillation  of  the  tank  makes  it  possible  to  determine  the  response  to 
an  arbitrary  forcing  function.  The  complete  solution  of  the  differential 
equations  of  motion  is  the  sum  of  the  complementary  function  and  the 
particular  integral.  The  complementary  function  corresponds  to  the  free 
oscillations  of  the  fluid  in  a stationary  tank,  and  the  desired  Initial 
conditions  can  theoretically  be  imposed  by  proper  choice  of  the  Infinitely 
many  arbitrary  constants  associated  with  the  fundamental  and  harmonics.  The 
particular  integral  for  an  arbitrary  forcing  function  could  be  obtained  by 
superposition  of  solutions  for  different  forcing  frequencies  through  a 
Fourier  Integral  approach,  for  example.  In  practice  the  Laplace  transform 
methods  seem  preferable. 


CONCLUSIONS 


In  many  cases  it  is  practicable  to  replace  the  fuel  by  spring- 
mass  systems  for  purposes  of  dynamic  analysis.  These  spring-mass 
systems  correctly  represent  the  fuel  even  for  angular  motions  of 
the  tank.  When  many  types  of  tank  motion  must  be  simultaneously 
considered  these  equivalent  systems  may  become  complicated,  and  pos- 
sibly lose  physical  significance  through  the  necessity  for  negative 
moments  of  inertia. 

It  seems  probable  that  masses  corresponding  to  higher  harmonics 
of  the  fuel  motion  may  often  be  neglected  (or  considered  as  fixed 
mass).  This  type  of  approximation  is,  of  course,  very  important 
in  simplifying  calculations,  and  deserves  further  study. 

When  equivalent  mechanical  systems  are  introduced,  any  desired 
method  may  be  used  for  solving  the  complete  dynamics  problem.  Laplace 
transform  methods  may  be  used  in  conjunction  with  the  mechanical  systems 
or  without  reference  to  them.  In  either  case  terms  corresponding  to 
higher  harmonics  can  probably  be  omitted. 
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DIAGRAM  OF  MECHANICAL  SYSTEM 
REPRESENTING  RESPONSE  OF  FUEL 
TO  HORIZONTAL  AND  PITCHING 
MOTIONS  OF  THE  TANK 


MASS  RATIOS  OF  THE  MECHANICAL  SYSTEM 
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RATIO  OF  FIXED  MASS  AND  FIRST  TWO  MOVING 
MASSES  TO  FUEL  MASS  vs  TANK  ASPECT  RATIO 
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FIGURE  6 
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NON-DIMENSIONAL  SPRING  CONSTANTS 
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FIGURE  10 
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